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^ ! Abstract 

■ Using higher-order derivative with respect to the parameter, we will give lower 
^ . bounds for variance of unbiased estimators in quantum estimation problems. This 
Q I is a quantum version of the Bhattacharyya inequality in the classical statistical 

' estimation. Because of non-commutativity of operator multiplication, we obtain 

three different types of lower bounds; Type S, Type R and Type L. If the parameter 
is a real number, the Type S bound is useful. If the parameter is complex, the 
Type R and L bounds are useful. As an application, we will consider estimation 
of polynomials of the complex amplitude of the quantum Gaussian state. For the 
fSJ , case where the amplitude lies in the real axis, a uniformly optimum estimator for 

I the square of the amplitude will be derived using the Type S bound. It will be 

shown that there is no unbiased estimator uniformly optimum as a polynomial of 
annihilation and/or creation operators for the cube of the amplitude. For the case 
where the amplitude does not necessarily lie in the real axis, uniformly optimum 
estimators for holomorphic, antiholomorphic and real-valued polynomials of the 
amplitude will be derived. Those estimators for the holomorphic and real- valued 
cases attains the Type R bound, and those for the antiholomorphic and real-valued 
' cases attains the Type L bound. This article clarifies what is the best method to 

^ ■ measure energy of laser. 

■ J. 1 Introduction 

■ Quantum estimation is an important theory in quantum information [7, 10]. It is not 
merely useful for many purposes, for example, evaluation of realized quantum informa- 
tion processing, but also is a fundamental problem in its own right. In this theory, we 
consider an optimization problem of measurements estimating unknown state, with re- 
spect to a risk function under appropriate restriction. A typical case, we adopt in this 
article, is minimization of variance under unbiasedness condition. What physicists call 
an observable is an unbiased estimator, and the variance is equal to the mean square 
error if the estimator is unbiased. 

It has been known that the quantum Cramer-Rao inequality gives a lower bound 
based on first-order derivative and the Schwartz's inequality [7, 10, 18]. Although there 
is only one Cramer- Rao inequality for classical statistical estimation [12], there are sev- 
eral different inequalities due to the non-commutativity of operator multiplication; SLD- 
type, RLD-type and LLD-type inequalities formulated by operators respectively called 
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Symmetric Logarithmic Derivative, Right Logarithmic Derivative and Left Logarithmic 
Derivative [11]. Yuen and Lax [18] showed that, for the quantum Gaussian state model, 
the homodyne and heterodyne measurement, respectively, uniformly attains the SLD- 
type bound for the one-parameter model and the RLD-type bound for the two-parameter 
model. Nagaoka [14] showed that the SLD bound can be locally attained for all one- 
parameter models. For asymptotic settings, there are many arguments on the quantum 
Cramer- Rao bound [1, 5]. 

For the classical estimation problems, Bhattacharyya [2] has improved the Cramer- 
Rao inequality extending the order of the derivative. For the classical Gaussian distri- 
bution model with unknown mean parameter 9, the uniformly optimum estimator for 
any polynomial g{6) is made by the Hermite polynomials, and it attains the classical 
Bhattacharyya bound. See [15, 16] for details in classical cases. In the quantum esti- 
mation theory, Brody and Hughston [3, 4] defined a Bhattacharyya-typc lower bound 
generalizing the SLD for pure states, and they analyzed asymptotic property. 

In this article, first, we will propose three quantum Bhattacharyya inequalities for 
mixed quantum states with one real or complex parameter. Generalization of the SLD 
gives the Type S lower bound for the real parameter case, and generalization of RLD and 
LLD gives the Type R and Type L bounds for the complex parameter case. Second, as 
the application, we will consider the quantum Gaussian state model where the amplitude 
parameter 9 is unknown. If 9 lies in the real axis, the uniformly optimum estimator for 
9"^ attaining the Type S bound is a self-adjoint observable given as a superposition of the 
number counting operator and the homodyne operator. This is realized by squeezing fol- 
lowed by the number counting. For 9^, there is no uniformly optimum unbiased estimator 
written as a polynomial of the creation and/or annihilation operators. If 9 lies in the 
complex plain, we will present uniformly optimum operators for polynomials g{9) oi9 and 
9 (the conjugate). If g{9) is holomorphic, i.e., dg/d9 = 0, then the optimum estimator is 
given by the heterodyne measurement and it attains the Type R bound. If g{9) is anti- 
holomorphic, i.e., dg/d9 = 0, then the optimum estimator is also given by the heterodyne 
measurement and it attains the Type L bound. If g{9) is real- valued, i.e., g{9) = g{9), the 
optimum estimator is given by some polynomials of the annihilation/creation operators 
and it attains both Type R and Type L bounds. 

In quantum optics, the Gaussian state is a simple model of laser whose complex 
amplitude is fluctuated around 9 E C, and the energy of the laser is proportional to 
l^^p. Using the Type R and L inequality, we will see that the counting measurement is 
optimum. If 6' G M is previously known, the Type S inequality shows that a counting 
measurement after a squeezing operation is optimum. 

This article is constructed as follows. In Section 2, our problem will be formulated and 
a known proposition of the quantum Cramer-Rao inequality will be shown. In Section 
3, new theoretical results of quantum Bhattacharyya inequality will be presented. In 
Section 4, our theory will be applied to the quantum Gaussian state model. Appendix A 
is the proof for Section 3, and Appendix B is that for Section 4. 

The author thanks the referees for useful comments. 
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2 Setup 



Suppose that there is a quantum system with an unknown state. Consider a set of 
candidate states {pg} for the system parameterized by 6* e 0. In this article, it is 
assumed that these density operators are invertible and that © = R or 9 = C. We will 
use C as the parameter instead of 9 when we need to remark that the parameter may not 
be a real number. Our interest is to estimate the true value of g{9), where : © — > © is a 
smooth function. Since the probabilistic error is inevitable, we consider an optimization 
problem of estimation under some restriction. 

An estimator M for g{6) is a Positive Operator Valued Measure (POVM) taking mea- 
surement outcomes in ©. (In a strict definition, an estimator should take measurement 
outcomes in g{Q) D ©, but, for theoretical convenience, we adopt the weaker definition 
in this article.) 

The expectation (=average=mean) of M is 



If E[M] = g{9) for any ^ e 0, M is said to be unbiased. We adopt the variance of M as 
the risk function; the variance is defined, in this article, as 



An unbiased estimator with the minimum variance among all unbiased estimators at a 
point ^ G is said to be locally optimum at 6. If an unbiased estimator is optimum at 
any 6* G 0, it is said to be uniformly optimum. 

The lower bound for the variance of unbiased estimators has been given by using the 
Schwartz's inequahty and the first-order derivative with respect to 9. This bound is called 
the quantum Cramer-Rao inequahty. See [7, 8, 9, 10, 11, 18] for the proof and related 
topics. 

Proposition Assume that = IR, the variance of any unbiased estimator M for g{9) 
satisfies 



IfT := g'{9){J'^)^^L^ + g{9) (A scaler x is identified with x times identity.) is free of the 
parameter, then the Projection Valued Measure (PVM) taking measurement outcomes in 
M given by the self-adjoint operator T is the uniformly optimum unbiased estimator for 
g{9) and the equality holds for (1). 

Assume that © = C, the variance of any unbiased estimator M for g{Q {C, G ©) 
satisfies 



IfT:= g'{C){J^) ^L^~^giC) /^ee of the parameter and ifT is normal, i.e., TT^ = T'^T , 
then the PVM taking measurement outcomes in C given by the spectrum decomposition 






(1) 



v[M] > WiOWJ"", 
v[M] > WiOWJ"^. 



(2) 

(3) 
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ofT is the uniformly optimum unbiased estimator for g{() and the equality holds for (2). 
Similarly, ifT :— g'{Q{J^)'~^L^ + g{Q is free of the parameter and is normal, it is the 
uniformly optimum unbiased estimator and the equality holds for (3). 

Here, J"^, L"^, J^, L^, and are defined as follows. 
Definition of and . For the case = R, let be a self-adjoint operator 
satisfying 

d peL^ + L^pe 

TeP' ^ 2 ' 

and then define as Ti[pe{L^Y]. is called Symmetric Logarithmic Derivative (SLD) 
and is called SLD Fisher information. 

Definition of L^, L^, and J^. For the case = C, let and be operators 
satisfying 

ift-Pfi". |/'< = iVc (5) 

where C := x+y/-[y, d/dC := {d/ dx-\f-id/ dy) /2 and d/dC, := {d/dx+\/^d/dy)/2 for 
real variables x and y. Then and are defined as Ti[pQL^[L^y] and Tr[L^p^(L^)''']. 

is called Right Logarithmic Derivative (RLD) and is called RLD Fisher informa- 
tion. Similarly, is called Left Logarithmic Derivative (LLD) and is LLD Fisher 
information. 

The definitions of L'^, and are not unique because the equations (4) and (5) have 
many solutions on a space where pg does not depend on 9. For example, the heterodyne 
measurement for the Gaussian model is obtained in the form [J^)~^L^ + 9 [9 & £) 
where is a solution in an extended system with an ancilla state, while the homodyne 
measurement {J^)~^L^ -|- ^ G M) needs no extension. In spite of the ambiguity of L', 
the inner product J' is uniquely determined. 

When the sample size is finite, there is no unbiased estimators uniformly optimum 
except for a few cases where the parameter space is flat with respect to the metric 
defined by [13]. When pe is not smooth with respect to 9, the difference instead of 
the derivative is useful [17]. 

There may be cases where the unbiasedness condition is so strict that no estimator is 
unbiased, and where the variance (or the means square error) is not appropriate as the 
risk geometrically. However, it is worth studying such problems with a view to gaining 
theoretical insight. 



3 Quantum Bhattacharyya inequality 

The quantum Cramer-Rao inequality is generalized by using higher-order derivative in- 
stead of the first-order derivative. 
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3.1 Quantum Bhattacharyya inequality of Type S for the real 
parameter case 

Consider the case = R. Let Lf := *(Li,L2, ...,Lk) be a column vector of self-adjoint 
operators satisfying 

peLk + Lkpe 

:Pe = ^ ■ (6) 



To simplify notations, we introduce a column vector 
differential operators, and we write (6) as 



\d/de,...,d}'/de^) of 



Dk[pe 



PuLl + Llpe 



Let Jjf be a A; X A; matrix where the (i, j')-th entry is 

Ji,j ■= ^^\peLiLj\. 

The definition of Jf is also simphfied as = Tr^xfefpe-^f *(-^f )] where Trj„xn[^] means 
taking the trace of each entry of an m x n matrix A, namely. 



TV. 



m.l 



/J 



/Tr[Ai,i] ••• Tr[Ai,„]\ 



\^Tr[A 



TO,1 



Though the definition of Lf is not unique for a system extension with a known ancilla 
state, Jf is uniquely determined. 

Assume that is invertible. The Quantum Bhattacharyya inequality of 'Type S' is 
given as follows. 



Theorem 1 If M is an unbiased estimator for g (9) , it holds that 

V[M] > \D,[g{9)]){J^)-'D,[g{9)]. 



(7) 



Especially, if T := ^Dk[g{9)]{J^y^L'^ + g{9) is free of the parameter, then the PVM M 
given by the self-adjoint observable T is the uniformly optimum unbiased estimator and 
the equality holds for (7). 

See Appendix A for the proof. 



3.2 Quantum Bhattacharyya inequality of Type R and Type L 
for the complex parameter case 



Consider the case = C. Let D'^ be a column vector 



d d d,^ rf2 ^2 



d'^ 



where 



d" 



dc-dc 



dC dC dC' dCdC dC''"' dC''''"' dC'-^dC' 



^ ( ^ J — 1 d\l^ f d / — ^ d 
2'^\dx dy) \dx^ dy 



d^ 
dC^ 
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The number of the entries is K :— k{k + 3)/2. Define column vectors and of K 
operators as the solutions to the equations 

Let and he K x K matrices given by 

:= TrK.K\pcL^{L^n := Tr^xi.[L^Pc(^^)1, 

where 

\Am,l ■■■ ^mj V4,n ^ij 

Applying this notation to D^, we define 

d d d? d? d'^ d^ 



ct _^ 

^ ■ \dC' dC dC' dCdC' dC"' dC^' dC'dC^-i' dC^ 



The definitions of and are not unique due to the system extension with a known 
ancilla state, but those of are unique. 

The Quantum Bhattacharyya inequalities of Type R and Type L are given as follows. 



Theorem 2 If M is an unbiased estimator for g{Q, then it holds that 



V[M]>D^^[giO]{J^)-'D^[g{a 



V[M]>D^'[g{C)]{Jt)-'Dl:[g{C)] 



L\-i nCr 



(8) 
(9) 



Especially, ifT := Df^diOK^k) ^-^fc +fl'(C) is free of the parameter and ifT is normal, 
then the PVM given by the spectrum decomposition ofT is the uniformly optimum unbi- 
ased estimator and the equality holds for (8). Similarly, ifT := Dl:"^ [g{()]{J^)~^ + g{() 
is free of the parameter, it is the uniformly optimum unbiased estimator and the equality 
holds for (9). 

See Appendix A for the proof. 



4 Application to the quantum Gaussian model 

For a known constant N > and an unknown parameter ^ e ©, let 



nN 



exp 



N 



\a){a\d'^a. 



Here, d^a means dxdy where a — x-\-\/^^y, and jo;) is the coherent vector of the complex 
amplitude a :— x -\- \J —\y, i.e.. 



W :=exp(-— )5:-^e 

n=0 * 
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where {en}^=o is the orthonormal system. 

The quantum Gaussian model is a generahzation of the classical model of Gaussian 
distributions, where the probability density is given as 



27r 

Here, ^ G M is unknown. In the classical estimation problem of 9'', the A;-th Hermite poly- 
nomial T(x) := (—l)''e^ ^'^{d'^ /dx'')e~^ is the uniformly optimum unbiased estimator 
which attains the classical Bhattacharyya lower bound. 

For the quantum Gaussian model, we consider two models; the real Gaussian model 
= M and the complex Gaussian model = C. For the real Gaussian model, we 
consider two cases; g{6) = 9^ and g(9) = 9^. For g{9) = 9"^, the optimum estimator is 
given by a PVM with measurement outcomes in 0(= M ^ 9{&) = {x \ x > 0}). For 
g{9) — 9^, it will be shown that any unbiased estimator given by an observable as a 
polynomial of the creation/annihilation operators can not be uniformly optimum. We 
will identify a self-adjoint operator with the PVM. 

Theorem 3 Suppose that = R. 

U fl'(^) — ^h,en the unbiased estimator 

(2iV+l)2^" ^+ (2iV+l)2 (2iV+l)2 

uniformly attains the Type S lower bound, so T is uniformly optimum. Here, a is the 
annihilation operator satisfying ae^ = v^Cn-i (^''^'d aa^ — a'^aga = I (identity). 

U 9{(^) = ^'^7 "0 unbiased estimator of the polynomial form of the creation/ annihilation 
operators can be uniformly optimum. 

See Appendix B.l for the proof. 

The value g{9) — 9"^ may be measured by the counting measurement a^a— constant, 
or by the square of the homodyne measurement (a -|- a^)^— constant. This theorem says 
that the optimum estimator is a superposition of these two measurements. We also note 
that this optimum estimator is realized as the counting measurement 6^6— constant after 
the following squeezing operation on the system; 

a\ 1 [N + 1 N \ f a\ f b 



For the complex Gaussian model, we will consider these three cases for a polynomial 

9{0 ofCeC; 

Holomorphic case: dg{()/d( = 0, 
Antiholomorphic case: dg{C,)/dC, = 0, 
Real- valued case: g{() = g{C). 

In each case, an optimal unbiased estimator will be presented by a PVM taking measure- 
ment outcomes in C. 
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A PVM M taking outcomes in C will be described by a normal operator T, that is, 

I uM{du). 
Jwec 

Since E[M] = Tr[pT], it holds that 

V[M] = Tr[p0{T - E[T]){T - E[T])t] = Tr[pe(T - E[T])\T - E[T])]. 

See [6] for details on normality and sub normality of operators for quantum measurement. 

For the holomorphic and antiholomorphic cases, we need to extend the system to 
describe the normal operators. Let /C be an ancilla system spanned by {/nj^^o; ^^'^ 
b be the annihilation operator satisfying bfn — v^/n-i bb^ — b^b — 1. Namely, the 
original annihilation operator a means a ® /, and the new one b may represent I <^ a. 
The original state is extended to /o/o- 

Theorem 4 Suppose that = C and g{9) is a polynomial of 9 and 9. 
If g{9) is holomorphic, the unbiased estimator 

T^g{a + b^) 

uniformly attains the Type R lower bound, so it is uniformly optimum. 
If g{9) is antiholomorphic, the unbiased estimator 

T^-g{a) + b) 



uniformly attains the Type L lower bound so it is uniformly optimum, where g{z) := g{z). 
If g{9) is real-valued, the unbiased estimator 

min(m,n) , / \ \ 

^-^ ^-^ Vmm m,n —r/ 

m,n r=0 ^ V ' / / 

^ min(m,n)! / a y+™'^^(0''"-H^^^+niax(o,^-n) 



r! VA^+L 

uniformly attains both Type R and Type L lower bounds so it is uniformly optimum, where 
See Appendix B.2 for the proof. 

For the holomorphic and antiholomorphic cases, the optimum estimators are realized 
essentially by the heterodyne measurement, that is, the estimated values for g{a + 6^) 
and g{a^ + b) are both obtained by operating g{-) to the heterodyne outcome. Hence 
they can be simultaneously carried out. On the other hand, for the real-valued case, no 
ancilla system is used so that it can not be measured simultaneously with the holomor- 
phic/antiholomorphic cases. 

For a real- valued case g{9) = Re (6')^, the optimum estimator is of the form (a -|- 
a^)^/4— constant, i.e., the square of the homodyne measurement. This measurement 
does not commute with that for g{9) — 9"^ {9 E R) of Theorem 3. 



Appendices 

The proofs for Theorems 1 and 2 will be given in Appendix A, and those for Theorems 
3 and 4 will be given in Appendix B. 
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A Proofs of Theorems 1 and 2 

The first lemma implies that, for any POVM estimator for g{0), there is a PVM which 
has the same expectation and a smaller variance. 

Lemma 1 Assume that M is a POVM taking measurement outcomes in ©. Let 

I uoM{du). 



Then, it holds that 



f \uj\''Tr[peM{doj)] > Tr[peTT^], (11) 

I \u\^T^\peM{du)\ > Tr[peT^T]. (12) 
Ju>ee 



Proof. The first formula (11) is obtained by 

\u;\^TY[peM{du;)] - Tr[p,rrt] 

pe [ {uj -T){uj -T^)M{d 
'- J wee 



= Tr 

Similarly, (12) is obtained by 

^TT[peM{duj)] - Tr[p,TtT] 
I /" (a) - T^) {u - T)M{du) \ > 0, 



I 



\U!\ 



wee 



Tr 



P6 



□ 

Therefore, for the proofs of Theorems 1 and 2, it is sufficient to show that, for the 
case e = M, if Trfp^T] = g{9) holds for any ^ e O then 

Trb(r - g{e)r] = Tr[poT'] - g^Of > 'D,[gmJl)-'D,W)i 
and, for the case O = C, if Tr[p^T] = ^(C) and Tr[p^Tt] = g{C) hold for any C e © then 

Tr[pc(r-^?(C))(Tt-^?(C))] 

= Tr[p^TTt] - |^?(C)P > Dl\g{C,)]{J^)-'Dl[g{C,)l 

npdT'-gmT-gm 

= TV[p,rtr] - \g{0\' > D^\g{Om-'D^[g{0]. 

If = M, T is a self-adjoint operator, for which the existence of the POVM is trivial. 
On the other hand, if © = C, one should consider normality and/or subnormality of T 
with extension of the system. See [6] . 
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Proofs of Theorems 1 and 2. 

The Theorem 1 for = M and the Theorem 2 for = C are proved by using the 
Schwartz's inequahty in a similar way. 

If, for any ^ e R, an self-adjoint operator T satisfies Trfp^T] = g{0), then 

^.ApeLliT - gm = TV,,, [^£^i±^(T - g{e))\ 

Let and be column vectors of A; + 1 operators and k + 1 scalers, respectively, given 

as 

Since 

- Tr \n (U'\^ - ( ^[^^ 

T .-TVfe+ixfc+iip^C/ [U )\- ^^^[^(^)] js ) 

is non-negative where V[T] :— Tr[pgT^] — g{Oy, it holds that 

Hence we obtain the Theorem 1. 

If, for any ( E C, an operator T satisfies Tr[p^T] = g{C), then 

Tr^[pc(T - giCML^y] = Tr^[(T - ^7(C))Pc:(^^)1 
= Tr4Df[p,]T] - giOTvKiDfipS = D^\g{0] 

where K ^ k{k + 3)/2. Let and be column vectors of X + 1 = (A; + l)(fc + 2)/2 
operators, and let and W'" be column vectors of X -|- 1 scalers, given as 

Let Vi := Tr[p^TTt] - and V2 := Tr[p^TtT] - \g{C)\^. Since 

are non-negative, it holds that 

{Wy{T^)-'W^ ^V,[T] - Df[g{C)]{J^)-'D^[g{0] > 
and {WyiT^r'W'^ ^V,[T] - D^\giC)]{J^)-'D^[g{0] > 0. 

Hence (8) and (9) of Theorem 2 are satisfied, respectively. □ 
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B Proofs of Theorems 3 and 4 

Lemma 2 If P{a) := exp{—aa/N) for a e C, then 

(a-Cr{a-CrP{a-C)\a){a\ 
( N 



■ N 



\N + 1J ' 

Proof. Since a = Yl°lo V^^i-i^ 
(a I at = a{a\. Moreover, since 

da 2 

d , . a 



crp(a-c)i«)(«i 


(13) 








= a\a) and 






1=0 * 





we have 



at|Q;)P(Q;) = exp(— q;q;/2) >^ eiA.iP(a) 

=|„)(a-A)p(„), (15) 

(Q;|aP(Q;) = exp(— q;q;/2) el,iP(Q;) 

= |q;)(^q;- ^)p(q;) (16) 
(for any smooth function P{a)), Recursively using these rules (15) and (16) with 

^^Pio^ - C) = ^^(« - and ±P{a - () = ^^(« - 0- 
we obtain the results (13) and (14). □ 

Lemma 3 Let p and q be non-negative integers and let n := _|_ g _|_ x_ jj 

p < q, the n-th entry of a Type R operator is 

min(p,g) 



y- r_-^^nnn(J>,g)-r ( max(p, g) \ min(p,g)! 

ymin(p, q) — rj r\ 



(^d _ Qr+ma.x(0,q-p) (^^ij _ Qr+max{0,p-q) 



The n-th entry of a Type L operator is 

min{p,(j) 



ri\rain(p,g)-r f niax(p, g) \ mm{p,q)\ 
Vmin(p, q) - rJ r! 



(^d — Qr+inax(0,q-p) _ Qr+inax{0,p-q) 
^ Nl{N + l)^-+max(0,p-g) ^^^^ 
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Proof. li p < q, the higher order derivative {dF'^'^ / d^C,cPC,)pc_ is calculated as 

1 /• , , , , dP+'i ( {a-Q{a-C,y 



/ \a)ia\ ,^ - exp - ^\d^a 

1 [ , ,dp f a-C y ( (o-CKo-Oa ,2 



1 /■ , , , , /(i^; d\p /' z 
''ttN 

{z (a-C)(a-0/N) 



p 

X 

r=0 



Similarly, if p > g, 



By applying Lemma 2 to (19) and (20), the n-th entry of is obtained as 



Z^r=OV -^J \q-r) r\ NP(N+iy P ^ <d, 



which means (17). Likewise, the n-th entry of is 



l^r=0\ ^) \p-r) r\ Ni{N+lY n \ (/, 

( l)<^~r( P \q\ (a-<;na^-(r-''+^ H V > Q 
.Z^r=OV \q-r) r\ Ni{N+1)p-i+^ P ^ 



which means (18). 



Lemma 4 Suppose that p, q, r, s are non-negative integers and that p + q and r + s 
not larger than k. Let 

{p + q-l){p + q + 2) {r + s-l){r + s + 2) 
m := hg + 1, n := hr + 1. 

Then the {m,n)-th entry of is 

exp 



dKPdKidX'dX^ \N N + 1 



K=\=0 



P'-Q'- 



NP{N+i)i if Tn — n, 
ifm^n, 



12 



and the {m,n)-th entry of is 



^p+q+r+s ^ RX 

exp 



dnPdKidX'^dy 



N +1 N 



K=\=Q 



Ni{N+i)p if n, 
if m ^ n. 



(22) 



Proof. Since is invertible under the assumption N > 0, — is a Type R 

operator for any C E C We have 

and hence the (m, n)-th entry is 

fp>+q+r+s 



K=\=0 



d^dR^dX^dXr^^^^'^'^^''^'^^'^ 

The parameter ( in (23) can be set to zero because 

Tr[p-Vc+KPc+A] = Tr[Up,'U-'Up,U-'UpxU-'] = Tr[po V.Pa] 
where U :— exp(Ca^ — (a). This can be calculated as 



1/3 -AP 



a 



n=0 

|/3)(/3|exp( - 



[ [ exp(-^- 



i'^Ny J^^c J pec 
— aa — PP + al3)d^ad^ (3 
N + 1 



1/3 -Ap + 1 



N 



+ 



N 



-a(3 



OC /'OC' /-OC /'OC 



{TrNy 
kX 



■oo •/ — oo •/ — oo ^ — oo 

rx 



exp ( - *(^; - Q-'p)Q{v - Q'V) 



) dxdydzdw 



where V :— *{x,y, z,w), p := *(Re(/t), Im(/t), Re(A), Im(A))/iV and 



Q:- 



[ l/N + l -l/(2iV)-l -v^/(2iV)\ 

1/iV+l v^/(2iV) -l/(2Ar)-l 

-l/(2iV)-l x/^/(2A^) l/N + l 

V-a/^/(2A^) -1/(2A^)-1 l/N + l J 



Since Q > and det g = (AT + l)7iV^ (24) is equal to 



f kX kX 
exp — + 



\N N+l 



(23) 



(24) 
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so the result (21) is obtained. 

Similarly, can be given as — Df[p(^]p^^. Hence 

so (22) holds. □ 

B.l Proof of Theorem 3 
Optimality of (10) for the case g(0) = 0^ 

First, we formally extend the real parameter 6 of pe to the complex parameter ( = 
9 + -\/^77 e C. Then the derivative [d/dOYpe can be considered as {d/C, + d/dC,Yp,^. A 
Type S operator Lf = *(Li, L2) is given as a solution to the equation 

d^ (d d\^ peLk + LkPe 

df^'' ^[Tc + dc)''^ 2 ^^'^ 

and Lk — L\ for /c = 1, 2. Let 

Lk-=Y^ Ci,j {a'a^^ + a^a^') (qj e R) . 

Since each coefficient for a'a^"* in the equation (25) should be zero, the solution is obtained 



as 



_ + 1 a-e N gt - g 

~ 2ivTTivTT^ 2Ar + l N 

N + 1 N 

=2 Mn 1 + 2 Ml 0, 

2N+1 2iV + l 

^ _ 2{a - ef + {a^ - Of 2{a-e){a^-e) 2 

^~ Ar2 + (7V + l)2 + ]V(]VTT) TV 

2(A^ + 1)2 (a-eV^ 2Ar2 /at_^x2 



TpIat + J ^ Ar2 + f at + 1)2 (, TV J 



A^2 + (Ar + i)2VAr + i; 7v2 + (Ar + i)5 
. (a_g)(at_g) 1. 

V N{N + 1) n) 
2{N + lf 2A^2 



where Mjj are Type R operators out of = *(Mi,o, -^^2,0, -^1,1, -/W"0, 2). As the 
inner product J 2 of is given in Lemma 4, is obtained as 



J; 



s 



( 2N+1 ^ ^ 

I ^ + ^ / ■ 



Since T of (10) is equal to (26',2)(J|^) ^Lf + 6*^, it attains the equality (7), hence it is 
uniformly optimum. □ 
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Non-existence of uniformly optimum estimator for the case g{0) = 0^ 

Since, for all non-negative integers m,n, the leading term of Tr[p5ia"*(a^)"] is 6"^6^, the 
form of an unbiased estimator for 6^ of a polynomial form of the creation/annihilation 
operators is given in the form 

T = u{a^ + a'^) + v(a^^a + a^a^) + w{o? + a^^) + xa)a + y{a + a^) + z. 

Using the characteristic function 

Tr[pee^'^^e^"] =\ I exp (\a + Aa - ^^-^ 
-exp((A + A)^ + 7V|A|2), 

we have 

Tr[per] = 2e^u + (2^^ + 4^(A^ + 1))^; + 2e'^w + (^^ + + l)x + 2^^ + ^. 
The unbiasedness condition tbr[p0T] = 9^ requires that 
1 X 

--V, --, y = -2(A^ + 1)^, z = -{N + l)x. 
It will be shown that, for any fixed ^ e ©, the variance is minimized if 

u———— — -T, v — ii- 



2(4A^2 + + 3) ' 2(4A^2 + 47V + 3) ' 

36 

W — ———Z --Z TT, X 



2{2N + l)2(4Ar2 + 4Ar + 3) ' (2A^ + lf{AN^ + 4Ar + 3) ' 

^ (N^ + N+im + l) ^ ~3{N+1)9 
^ Am + 4N + 3 ' (2iV + l)2(4Ar2 + 4iV + 3) ■ ^ ' 

Since w, x and z depend on 9, there is no unbiased estimator uniformly minimizing the 
variance. 

By the same way as the previous proof, the third entry of Lf = *(Li, L2, 1/3) can be 
obtained as 

{a-9f + {a^-9f (a - 9y{a^ - 9) + {a - 9){a^ - 9^ 
^~ N^ + {N + iy ^ N{N + 1){2N + 1) 

a + a^ - 29 
N{2N + 1) 

2{N+lf {a-9f 2N^ {0) - 9f 



+ (AT + 1)3 (TV + 1)3 7V3 + (Ar + i)3 

Q{N + lfN / {a-9Y{a) -9) _ 2{a - 9) \ 
^ N{N + 1){2N + 1) \ {N + 1)^N {N + 1)n) 
6{N + 1)N^ ({a-9){a^~9f 2(at - 6 



+ 



iV(iV + 1)(27V+ 1) V (Ar+l)7V2 m 

2{N+lf 2N^ 

I 6(Ar + i)^Ar , 6(Ar+i)Ar2 
Ar(A^ + l)(2A^ + l) Ar(Ar + l)(2Ar + l)''^''^ 
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where Mjj are Type R operators out of = *(Mi,o, -^3,0, ^^2,1; ^1,2, M),3)- Using 
Lemma 4, we have 







= - 

n .,s,?t,..s + 



Ar3+(Ar+l)3 ^ Ar(Ar+l)(2Ar+l), 

The unbiased estimator satisfying formulas in (26) is equal to T = *D^[9^][J^)~^L^ + 0^, 
which attains the lower bound (8) for the variance at each 9. □ 

B.2 Proof of Theorem 4 
The holomorphic case dg/dd — 

Consider the monomial case g{9) — The column vector -D^ [C*^] is given as 

k\ 
J 

Lemma 4 shows that is a diagonal matrix of the form 



^^[C'] = * (0, K^-\ 0, 0, k{k - 0, 0, 0, 0, k\) . 



, / 1 1 2 1 2 3! 

-^k "'^^Vat' AT + l' Ar2'iV(iV+l)' (Ar + l)2'iV3 



' Ar'^ ' ■ ■ ■ ' 



(j-1)! j\ (j + 1)! (A:-l)! k\ 



N{N +1)0-^' {N + 7Vi+i A^(iV+l)'=-i' (AT + l) 



k y • 



The system is extended with the ancilla space /C, and the state is set as (8) /o/o- The 
annihilation operator on K, is denoted by h. A Type R operator on the extended 
system is 



V AT ' AT + i ' ' Ar(Ar + i) ' (at + 1)2 ' 

(gt - C)^ (a + bt_^)j-i(at_^) (g + &t - C)^- (gt - 0^+^ 



N{N + lY-^ ' (AT + l)*^ 
Define an operator as 

A; 



Since IV^xik^f] = ^fe[Tr[pc]] = = 0, Tr[p^Tfc] = (K Hence is uniformly 

optimum unbiased estimator for C*^. For the general holomorphic case g{() — J^^. Ck9'' 
(cjfc e C), the estimator T :— g{a + b^) — 'Ylik^kTk is unbiased and uniformly optimum. □ 
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The antiholomorphic case dg/dd = 

Consider the monomial case g{0) = The column vector [C^] is given as 

J 

By Lemma 4, Jj^ is a diagonal matrix of the form 



^."^[ci = * {ke-\ 0, k{k - l)c^-^ 0, 0, 0, o) . 



4 -diag(^^^^,^, (jv + i)2' Jv(Jv + l)'••■' 
(i-l)! j! (j - 1)! 

For the annihilation operator h on the ancilla system for the extension /o/o; a- TyP© 
L operator is 



V AT+i ' AT ' (iv + l)2 ' iV(iV+i) 

ATJ-i ' (iV+l)J ' iV(iV + l)i-i 
{a-Cf-\a^ + h-Q (a -OS 



iVfc-i(Ar + i) ' ATfc 
Define an operator as 

k 



j=o 

Since Tr^xi[pc^fc] = ^fe[Tr[Pc]]_= = 0, Tr[pcTfc] = C'- Hence n is uniformly 

optimum unbiased estimator for C^^. For the general antiholomorphic case g{Q = CkO^ 
(cjk e C), the estimator T := g{a^ + b) = ^j^. CfeTfe is unbiased and uniformly optimum. □ 
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The real-valued case g = g 

By Lemmas 3 and 4, T^,„ := Df[cm,nC"'C]iJk)~'L^+n + for A; > m + n is 



p=U q 



q \p!(a-C)^-^+'-(at-C)' 



r=0 



p-rj r\ Np{N + l)<}-p+r 



rjr\ NP{N+iy 



ip<q) 

{p > q) 



{m — q)\{n — p)\ 



X < 



p=0 q=0 

p 

r=0 

^ (-l)«-'-(a-C)''(a^ -0^""+" 

I, r=0 



{p-r)\{q-p + r)\r\{N + 



(g - r)! (p - g + r)! r!(A^ + 



{p > q). 
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Tm,n does not depend on C £ C because dTm,n/ dC is equal to 



p=U q= 

P 



X < 



E 

'r=0 
9 



(-l)P-^(a - C)'^~^+''(a^ - 0' 
(jo - r) ! (5 - p + r) ! r!(A^ + !)-?'+'■ 



{p < q) 



I, r=0 



(g - r)! (p - ? + r)! r!(A^ + 
- WZ^Z^ (m-g)!(n-p)! 

p=0 q=l ^ ^ ^' 



X < 



r=0 



(p - r)! (g - p + r - 1)! r!(iV + 
^ (g - r)! (p - g + r)\ (r - 1)!(A^ + 1)-"^+ 



{p<q) 

{p>q) 



n m— 1 



0^ (m-g- l)!(n-p)! 

(-l)f-'-(a - C)^-^+^(at - 0' 



X < 



P 

E 

r=0 
5 



{p-ry.{q-p + r)\r\(N + 1)-p+^ 



2-^ (n- 



( r=0 



{q-ry.{p-q + ry.r\{N + 



ip < q) 

(p > q). 



n m— 1 



m!n!C 



im — q' — l)\in — p)\ 



( P 



^ (-l)P-"(a-C)^'"^+"(a^ -0' 

Z-^ in — 



X < 



r=0 



(p - r)! (g' - p + r)! r!(A^ + 1)- 



p+r 



^(g'-r')!(p-g' + r')!r'!(iV +!)-«'+'■' 



=0, 
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where q' :- 
to 



q — 1 and r' :— r — 1. Tm,n does not depend on C, because dTm^n/dC, is equal 



li— 1 in 
p=0 q=Cl 



m!n!C™"''C 



m—q/^n—p—l 



{m ~ q)\{n — p — 1)\ 



X < 



E 

r=0 



{p -ry.{q-p + r)\ r\{N + 1)-p+^' 



(g - r)! (p - g + r)! r!(7V + l)-</+'- 
(m — — p)! 



(p < q) 



{p > q). 



EE 



p=l q=0 



E 



X < 



E 



^ (p - r)! (g - p + r - 1)! (r - 1)!(A^ + 1)-?'+'^ 



I, r=0 

n— 1 m 



(g - r)! (p - g + r - 1)! r!(A^ + 



(p < g) 
(p > g) 



E 



(m — g)!(n — p — 1)! 



X < 



(p - r)! (g - p + r)! r!(iV + 
^ (-l)^-^-(a - C)'-(at - C)P-''+'- 



n— 1 m 



(g - r)! (p - g + r)\ r\{N + 1)" 



(p < g) 



X < 



=0, 



^(p'-r')!(g-p' + r')!r'!(iV+l)-P'+'^' 



E 

r=0 



(g - r)! (p' - g + r)! r\{N + 



(p' > g) 
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where p' 



p — 1 and r' :—r — 1. Therefore, 



T 



m,n 



\ I — " 

min{m,n) 

=Cm,n{N+ir E (-1) 
r=0 




{q - ry. ip - q + ry. r\{N + 



;^ (p-r)!(g-p + r)!r!(iV+l)-P+'' 



{-l)P-''a'i-P+''{a^Y 




{p < q) 




min(m,n)!/ a \ r-+max(0,n-m) 

I r+max(0,m— n) 




Let g{C) = Xlr„,„ be a real- valued polynomial, namely, Cm,n = c^. Let T := Y.m,n Cm,nTm,n% 
Then, T is a self-adjoint observable. Since Trfexi[pc-^fe] — 0' Tr[p^T^^„] = c^.tiC^C"- 
Hence T is an unbiased estimator for g{0- T attains the Type R lower bound because 



It can be shown in the same way that T = D^[g{Q)]{Jj^)-^L]^, so that T attains the 



References 

[1] Asymptotic Theory Of Quantum Statistical Inference: Selected Papers (Edited by 
Masahito Hayashi), World Scientific, New Jersey. 

[2] Bhattacharyya, A. 1948 Sankhyd 8, 315-328. 

[3] Brody, D. C. and Hughston, L. P. 1998 The geometric universe (Oxford, 1996) 265- 
276, Oxford Univ. Press, Oxford. 

[4] Brody, D. C.and Hughston, L. P. 1998 R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. 
Sci, 454, 2445-2475. 

[5] Gill, R. and Massar, S. 2000 Phys. Rev. A 61, 042312. 

[6] Hayashi, M and Sakaguchi, F. 2000 J. Phys. A 33, 7793-7820. 

[7] Helstrom, C. W. 1976 Quantum Detection and Estimation Theory Academic Press, 
New York. 

[8] Holevo, A. S., 1973 Journal of Multivariate Analysis 3, 337-394. 
[9] Holevo, A. S., 1973 Teor. Verojatnost. i Primenen., 18, 371-375. 
[10] Holevo, A. S. 1982 Probabilistic and Statistical Aspect of Quantum Theory, North- 



[11] Holevo, A. S. 2001 Statistical structure of quantum theory. Lecture Notes in Physics. 
Monographs 67. Springer, Berlin. 



Type L lower bound too. 



□ 



Holland. 



21 



[12] Lehmann, E. L. 1983 Theory of point estimation Wiley, New York. 
[13] Matsumoto, K. 2002 J. Phys. A 35, 3111-3123. 

[14] Nagaoka, H. 2005 Asymptotic Theory Of Quantum Statistical Inference: Selected 
Papers (Edited by M. Hayashi), 113-124, World Scientific, New Jersey. 

[15] Tanaka, H. 2003 Comm. Statist. Theory Methods 32, 1885-1896. 

[16] Tanaka, H. and Akahira, M. 2003 Ann. Inst. Statist. Math., 55, 309-317. 

[17] MR21 17272 (2005j:81026) Tsuda, Y. and Matsumoto, K. 2005 J. Phys. A 38, 1593- 
1613. 81P15 (94A15) 

[18] Yuen, H. P., and Lax, M. 1973 IEEE Trans. Information Theory IT-19, 740-750. 



22 



